Multi-spin dynamics of the solid-state NMR Free Induction Decay 
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We present a new experimental investigation of the NMR free induction decay (FID) in a lattice 
of spin-1/2 nuclei in a strong Zeeman field. Following a n/2 pulse, evolution under the secular 
dipolar Hamiltonian preserves coherence number in the Zeeman eigenbasis, but changes the number 
of correlated spins in the state. The observed signal is seen to decay as single-spin, single-quantum 
coherences evolve into multiple-spin coherences under the action of the dipolar Hamiltonian. In 
order to probe the multiple-spin dynamics during the FID, we measured the growth of coherence 
orders in a basis other than the usual Zeeman eigenbasis. This measurement provides the first 
direct experimental observation of the growth of coherent multiple-spin correlations during the FID. 
Experiments were performed with a cubic lattice of spins ( 19 F in calcium fluoride) and a linear 
spin chain ( 19 F in fiuorapatite). It is seen that the geometrical arrangement of the spins plays a 
significant role in the development of higher order correlations. The results are discussed in light of 
existing theoretical models. 
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Solid-state NMR is an ideal test-bed for studying co- 
herent quantum dynamics in a large Hilbert space. In 
this work, we experimentally investigate the many-spin 
dynamics of the NMR free induction decay. The free in- 
duction decay (FID) is the response of the spin system 
following a tt/2 pulse. In a solid lattice of spin-1/2 nuclei 
in a strong magnetic field, this evolution is dominated by 
the secular dipolar Hamiltonian. This is a classic prob- 
lem in spin dynamics and has been well-studied since the 
early days of NMR 1-11 . At the magnetic fields typically 
used in NMR experiments (< 20 T), the spin system is 
highly mixed at room temperature, and its equilibrium 
state is represented by a thermal density matrix. If the 
external Zeeman field is much stronger than the inter- 
nal dipolar fields of the sample, the normalized deviation 
density matrix in thermal equilibrium can be approxi- 
mated as 



p(o) = -E4 



(i) 



Following a 7r/2 pulse, the spins evolve under the secular 
dipolar Hamiltonian 
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The strength of the dipolar coupling Djk between spins 
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j and k is given by 



D jk = 
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where 7 is the gyromagnetic ratio, Tjk is the distance 
between spins j and k, and 6jk is the angle between the 
external magnetic field and internuclear vector fjk- Since 
the Hamiltonian is time- independent, the formal solution 
of the Liouville-von Neumann equation yields the density 
matrix of the spin system at time t following the pulse as 
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An exact solution to this many-body problem has not 
been found, but the equation can be expanded in a power 
series to examine the short time behavior of the system: 



/5(0) + -t p{Q),H mt 



2h 2 



• Hu 



(5) 



In an inductively detected NMR experiment (in which a 
coil is used to measure the average magnetization), the 
observed signal is given by 



S(t)=C{U)=CTr{Up(t)} 



(6) 



where 7+ = ^2j(Ij X + Hjy) an d C is a proportionality 
constant. The only terms in pit) that yield a non-zero 
trace in the above equation and therefore contribute to 
the observed signal S(t) are the single-spin angular mo- 
mentum operators such as which are single-spin, 
single-quantum coherences. Single quantum coherences 
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are off-diagonal terms of the density matrix (in the Zee- 
man eigenbasis or the z-basis) connecting eigenstates 
with Am = ±1 (corresponding to coherent superposi- 
tions of these eigenstates). Evaluating the commutators 
in Eq. (5) 



P(t) 
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Substituting Eq. (7) into Eq. (6), it can be seen that 
the observable magnetization decays during the evolu- 
tion under H- mt because single-spin, single-quantum co- 
herence terms are transformed to unobservable multiple- 
spin, single-quantum coherence terms by the higher-order 
nested commutators. The n th term in the expansion in 
Eq. (5) has n-spin correlations. 

There has been much theoretical effort to predict the 
shape of the FID in calcium fluoride (CaF 2 ) 1 ~ 8 . Calcium 
fluoride is a standard test system for spin dynamics as the 
19 F (spin-1/2) nuclei are 100% abundant and form a sim- 
ple cubic lattice. The main goal has been to reproduce 
the decay and beat pattern of the observed time domain 
NMR signal. For example, Engelsberg and Lowe 10 mea- 
sured up to 14 moments of the FID in CaF2, and these 
were found to be in good agreement with theoretically 
calculated values for the 2nd to 8th moments. The odd 
moments of the FID are zero, and the even moments are 
given by 1 
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Evaluating the nested commutators becomes increasingly 
challenging and the higher-order moments are difficult 
to calculate. However, it is these higher moments that 
characterize the many spin correlations in the spin sys- 
tem. It can be seen that the 2n th moment arises from 
the (2n+l) th term in the expansion in Eq. (5), which 
creates up to (2n+l) correlated spins. The main weak- 
ness of the moment method lies in the fact that the most 
important contribution to the value of the higher mo- 
ments comes from the tails of the FID, which are acquired 
with the lowest signal-to-noise ratio (SNR) in typical FID 
measurements 3 . 

In this paper, we use a novel multiple-quantum NMR 
technique 14 to study multiple-spin dynamics during the 
FID. Standard multiple-quantum techniques 15-18 encode 
coherence orders in the Zeeman eigenbasis (or z-basis), 
but coherence numbers are conserved under the secular 
dipolar Hamiltonian in this basis 19 . In our experiment 
we encode multiple-quantum coherences in the x-basis. 



The dipolar Hamiltonian in the x-basis is 
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and no longer conserves coherence order in this basis (we 
use the superscript x to denote that the raising and lower- 
ing operators are defined in the x-basis; these operators 
are otherwise assumed to be expressed in the z-basis). 
The coherence orders are encoded by a collective rotation 
about the x-axis (which is the effective quantizing axis 
in this basis). Transforming the density matrix shown in 
Eq. (7) into the cc-basis yields 



p x (t) = -E f i--! ft E D «(^-M-) + 
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From Eq. (10), it can be seen that even order multi- 
ple quantum coherences are created in the x-basis. It is 
possible to generate only odd order coherences using a 
y-basis encoding for the same initial state. 

It is useful to consider the dipolar evolution of this 
highly mixed state using the Liouville space formulation 
for multiple-quantum dynamics suggested previously 17 . 
The density operator in Liouville space can represented 
as 



N 



K 
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(11) 



K=0n=-K V 



where Pku P represents a basis operator which is a prod- 
uct of K single-spin angular momentum operators, n is 
the coherence order of the operator and p is a label that 
identifies a particular configuration of spins having the 
same K and n. The selection rules for the dipolar Hamil- 
tonian in the Zeeman basis are given by 



AK = ±1 , An = 



(12) 



A projection of Liouville space onto the two-dimensional 
plane spanned by K and n is shown in Figure 1(a). Fol- 
lowing a 7r/2 pulse, the trajectory in the Zeeman basis 
is indicated by the arrows (only positive coherences are 
shown here; the evolution is perfectly symmetric for neg- 
ative n). Increasing numbers of spins are correlated fol- 
lowing evolution under the dipolar Hamiltonian, but the 
coherence number does not change. Figure 1(b) shows 
the same evolution in the x-basis where the selection rules 



are 



AK = ±l, An = 0,±2 



(13) 
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(a) z-basis ( b ) x-basis 
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FIG. 1: Projection of Liouville space onto the two-dimensional plane spanned by K and n, showing the dynamics of the FID 
in (a) the Zeeman eigenbasis, and (b) the x-basis. The arrows show the allowed paths in each case. 
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FIG. 2: Pulse sequence used in this experiment, t is the 
evolution time under the secular dipolar Hamiltonian. The 
48-pulse sequence was used to suppress the evolution of the 
internal Hamiltonian during the i-basis encoding step. A 
magic-echo sequence was used to reverse the dipolar evolution. 
A delay (5=10 fjs) was inserted before the magic echo to push 
the echo out, in order to minimize any pulse transients and 
receiver dead time effects. The duration of the spin locking 
time in the magic echo sequence is 2ti oc i t where ii oc k = t + 
25 + 3/xs. 



Starting from an initial I x state (K = l,n = 0), only 
even order coherences are observed. In this paper we 
characterize the growth of these coherences. 

The pulse sequence used in this experiment is shown 
in Figure 2. After an initial tt/2 pulse, multiple-spin, 
single-quantum states in the Zeeman basis are created 
during evolution under the secular dipolar Hamiltonian 
as described in Eq. (7). A <fil x rotation encodes coher- 
ence orders in the x basis, and a magic-echo sequence 1 is 
used to refocus the multiple-spin terms back to observ- 
able single-spin, single-quantum coherence terms. The 
(j)I x rotation is obtained by applying two tt/2 pulses, 
with phases y + <fi and y, which results in the prop- 
agator exp(i<fil z ) exp(i<pl x ). The initial tt/2 excitation 
pulse is also phase shifted by cj) to cancel out the rota- 
tion about I z . Since it is difficult to apply back to back 
tt/2 pulses without a delay between them (without intro- 
ducing phase transients or allowing some dipolar evolu- 
tion during the pulses), an evolution suspension sequence 
needs to be used in between the two tt/2 pulses. In this 
experiment we use a previously described 48-pulse evo- 
lution suspension sequence 13 . 

The experiments were performed at room temperature 
at 2.35 T (94.2MHz, 19 F), using a Bruker Avance spec- 
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FIG. 3: i-basis coherence order distribution at various time 
points under the evolution of the secular dipolar Hamilto- 
nian in CaF2. The multiple-quantum intensities have been 
re-normalized to set the intensity of the zero-quantum term 
to one in each case. 



trometer and home-built probe. The samples used were a 
1 mm 3 single crystal of CaF2 with T\ ~ 7 s, and a crystal 
of fluorapatite (FAp) with T Y ~ 200 ms. The FAp crys- 
tal is a mineral crystal specimen from Durango, Mexico. 
Such natural crystals are expected to have a number of 
interruptions in the spin chains as well as a variety of 
other defects 20 as indicated by the short Ti. All experi- 
ments were conducted on resonance. High power 0.5 [1% 
tt/2 pulses were used for the 48-pulse evolution suspen- 
sion sequence, while lower power 1.5 /is 7r/2 were used 
during the magic-echo sequence as this sequence is more 
susceptible to phase transient errors. The phase <fi was 
incremented from to 4-7T with A(j) = tt/32 to encode up 
to 32 quantum coherences for every experiment. A fixed 
time point corresponding to the maximum intensity sig- 
nal was sampled for each <f> value and Fourier transformed 
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with respect to 4> to obtain the coherence order distribu- 
tion at each dipolar evolution time t. 

Figure 3 shows the coherence order distribution ob- 
served for CaF2 at various time points during the FID. 
At short times, the maximum coherence order (n max ) 
corresponds to the maximum number of correlated spins 
i-e Qk n P = for K > K max ). At longer times, 
the maximum coherence order observed in the experi- 
ment sets the lower limit of the size of the spin correla- 
tion, since the SNR of higher order coherences might be 
too low to be observed. 

Figure 4 shows the growth of the different coherence or- 
ders in CaF2 during the FID. The inset shows the initial 
oscillation between the zero and double quantum coher- 
ences at short times (which corresponds to single and two 
spin correlations respectively) due to the resolved near- 
est neighbor coupling at the [100] and [110] directions. 
This oscillation may be theoretically understood by con- 
sidering the time development of an isolated pair of spins 
under the secular dipolar Hamiltonian (in the x-basis): 
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where D is the strength of the pairwise coupling. In an 
extended spin system, this oscillation is rapidly damped 
by leakage from isolated pairs to higher order correla- 
tions. The higher order coherences (n >4) are seen to 
follow a sigmoidal growth curve. The higher order coher- 
ences develop later in time, and this progressive growth 
leads to a saturation of the intensities of the lower order 
coherences, consistent with the model in Figure 1. As 
t increases, imperfect refocusing of the dipolar evolution 
under the magic echo sequence results in a decay of the 
observed signal. In order to remove this decay the inten- 
sity for each coherence order is normalized with respect 
to the total signal measured at that evolution time. 

We have fit the initial growth of each coherence order 
to the following sigmoidal function 



S n (t) = 



1 



,—a n (t—t° 



(15) 



In Figure 5 we plot the variation of a n and C n as a 
function of coherence order n. The parameter a n rep- 
resents the underlying rate at which the different coher- 
ence orders are transformed, and should be dominated by 
the strength of the dipolar couplings involved. It is seen 
that a n does not vary with n, suggesting that the near- 
neighbor interactions dominate the dynamics here. The 
mean dipolar coupling strength depends on the crystal 
orientation, so we measured a n with the crystal oriented 
along the [111], [110] and [100] directions with respect 
to the external field. The mean dipolar coupling can be 
estimated by summing \Du\ over the 26 nearest neigh- 
bor spins for each crystal orientation. For a simple cubic 
lattice, the ratio of these means is 1:1.45:1.87 for [111], 
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FIG. 4: The growth of multiple-spin correlations during the 
FID, showing the sigmoidal fit to the initial growth data of 
each re-basis coherence order (n > 4), along ~[110] direction in 
CaF2. The inset shows the dynamics of n — and 2. The in- 
tensity of the signal for each coherence order S n is normalized 
with respect to the total signal Stot, in order to compensate 
for imperfect refocusing of the dipolar interaction. 



[110] and [100] in good agreement with the ratios of a n 
shown in the figure. The values of C n are seen to decrease 
as coherence number increases, independently of crystal 
orientation. This is expected as the total polarization 
is conserved, and the signal is progressively spread over 
increasingly larger regions of the system Hilbert space. 

Figure 6 shows the onset time i° nset of each of the x- 
basis coherences for different orientations of the crystal. 
Physically, the onset time corresponds to the time re- 
quired for a specific coherence order to become observable 
in the experiment. While a first glance at Eq. (7) would 
seem to suggest that higher-order correlations should de- 
velop as i", it is the geometry of the spin system (the 
values of -Dy) which ultimately determines the rate at 
which the spin correlations grow. The onset times de- 
pend on the rate at which the correlations are spreading 
through the spin system, which in turn depends on the 
value of the dipolar couplings. Thus the rate is expected 
to be fastest (shorter onset time) with the crystal ori- 
ented along the [100] direction and slowest (longer onset 
time) for the crystal oriented along the [111] direction, 
in agreement with the experimental data. The variation 
of onset times with coherence order is also observed to 
depend on the dimensionality of the spin system. The 
variation is sub-linear in the cubic CaF2 system and dis- 
plays an approximate n 2 / 3 dependence. The inset in Fig- 
ure 6 shows that the onset times obtained for FAp, the 
quasi-one-dimensional spin system, vary linearly with co- 
herence number, in marked contrast to the results from 
CaF 2 . 

A variety of models have been proposed to describe 
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FIG. 5: The values of a n for different orientations in 
CaF2. The left inset shows the ratios a[ioo]/ a [m] (*) an d 
Q!~[iio]/ Q [iii] (+) • The calculated ratios of the mean dipo- 
lar coupling, obtained by averaging over 26 nearest neighbors, 
[100] /[111] (solid line) and ~[110]/ [111] (dotted line) are also 
shown. The right inset show the values of CVs for different 
orientation in CaF2 



FIG. 6: Onset time of multiple-spin correlations for different 
crystal orientations in CaF2 . The inset shows the onset times 
for the FAp sample. Odd coherence orders for the FAp sample 
were obtained using a y-basis encoding on the same initial 
state. The continuous lines represent the best fits of Eq. (20) 
to the data, assuming that the same equation is valid for 
coherence number as well. 



the dynamics under a multiple-quantum Hamiltonian. 
The most commonly used model involves a random walk 
among the components of the Liouville space basis set 
PKnp, subject to the selection rules of the multiple quan- 
tum Hamiltonian 17 . The model replaces the Liouville- 
von Neumann equation by a set of coupled rate equations 
with exponential solutions 

j f g = R-g , (16) 

where the vector g containes the coefficents gunp- All 
possible configurations that contribute to a particular co- 
herence are assumed to be present in equal measure, and 
the resulting growth of the spin system is described by 
a hopping procedure between the allowed points on the 
lattice (shown in Figure 1). Under this assumption, the 
hopping rates are solely determined by the degeneracies 
of the coupled states. The model thus eliminates any os- 
cillatory solutions and precludes the possibility of quan- 
tum interference effects playing a role in the evolution. 
The model also ignores the geometrical arrangement of 
the spin system and the specific distribution of dipolar 
coupling coefficients responsible for driving the evolution. 
All spin systems display a universal growth kinetics, as 
long as the dynamics are scaled by a lattice parameter 
that represents the mean dipolar coupling strength of the 
system, similar to the parameter a n obtained above. A 
numerical solution of the coupled equations for the mul- 
tiple quantum evolution was observed to yield sigmoidal 
growth curves for the various coherence orders 17 . Limi- 
tations to this model have been discussed elsewhere 21 . 



Munowitz and Mehring 11 used this model to numeri- 
cally simulate the growth of multi-spin dynamics of the 
FID in a 21-spin system. In order to track the devel- 
opment of correlations among the spins, they defined an 
induction time tx over which a particular X-spin coher- 
ence reaches half its maximum value. This parameter is 
very similar to the experimentally measured onset times 
of the different coherence orders described here. Figure 7 
of Ref. 1 1 shows the variation of induction time with the 
number of correlated spins. The numerical results show 
that the variation of the induction time is sub-linear for 
small numbers of correlated spins (< 10), in agreement 
with the experimental data. For larger numbers of cor- 
related spins, however, there is a marked deviation from 
the sub-linear behavior, as the effects of the finite system 
size (21 spins) begin to manifest themselves in the simula- 
tions. The number of correlated spins would have to ap- 
proach the number of spins in the sample (« 10 21 ) before 
such effects would be observed experimentally. While 
providing some insight into the growth of the dynamics 
for the cubic spin systems, it is seen that the model does 
not do a very good job at describing the one-dimensional 
spin system. It is in this situation that the spin geometry 
plays a dominant role. 

Gleason and co-workers have proposed an alternate 
model to describe the growth of spin correlations that 
emphasizes the geometrical ordering of spins 22 . By ag- 
gregating over the different configurations and coherence 
orders, i.e. summing over n and p in Eq. (11) the den- 
sity operator p is expressed as a sum of terms with spin 
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number K and coefficients gx, 



N 



P(t) = J2 9K(t)P K ■ 



(17) 



K=0 



where qkPk = J2„ P 9KnpPKn P - Essentially this model 
assumes a single effective if-spin operator that incorpo- 
rates all the possible spin and spatial configurations of 
the K spins. The resulting model for spin propagation 
through a lattice yields a differential equation for the co- 
efficients qk (t) of the form 



dt 



9k = -i {W K _ l9K -i + W K+1 g K +i) (18) 



where the rate constants W { and W T correspond to the 
forward and reverse rates respectively. Under the as- 
sumption that the spatial grouping of the K spins is con- 
tinuous, and that only the nearest neighbor couplings are 
important, the forward rate (and equivalently the reverse 
rate) can be expressed as W K oc Dn n n S} where D is the 
strength of the nearest neighbor coupling, n s is the num- 
ber of spins on the surface of the spatial grouping, and n n 
is the number of neighboring spins coupled to each spin. 
New spins are added on the surface of the correlated spin 
cluster. While n n is a constant, the term n s would differ 
significantly for spin systems of different dimensionali- 
ties, and can be expressed as n s oc K 1 ~ 1 / d where d is the 
dimensionality of the spin system. For a linear spin chain 
d = 1 and n s is independent of K, while for a cubic spin 
system, d = 3 and n s oc K 2 ^ 3 . While this model does 
not discuss coherence order, the dimensional dependence 
does agree with the experimental results, if the onset time 
characterizes the effective rate constant. In the limit of 
large K, the rate constants W t K _ 1 « W K+1 , and the co- 
efficients gait) are approximately given by 



g K (t) oc i 



K-l 



tanh 



(ptK- 1,d ) 



K 



(19) 



where /3 is proportional to the mean dipolar coupling. 
With appropriate choice of normalization, the magni- 



tude of gxit) as given by Eq. (19) strongly resem- 
bles the sigmoidal growth characteristics of the multiple 
quantum coherence intensities shown in Figure 4. An 
onset time can be obtained from Eq. (19) by setting 
[tanh{PtK- 1 l d )} K = 1/2, yielding 



t 1/2 = |i^/Vctanh (2-V*) 



(20) 



where r\ is a constant scaling factor. Figure 6 shows the 
best fit of Eq. (20) to the experimental data, assuming 
that the model holds true for coherence number as well. 
It is seen that there is excellent agreement at larger val- 
ues of n for the cubic CaF 2 system. The values of rj/P 
obtained from the fit are 37.47 in the [111] direction, 
31.13 in the [110] direction and 21.22 in the [001] di- 
rection. Their inverses are in the ratio 1:1.46:1.76 for 
[1 1 1] : [1 10] : [100] which is in excellent agreement with the 
theoretically calculated values shown earlier. For the lin- 
ear FAp system, Eq. (20), which is linear for large K, is 
observed to be weakly superlinear at these small values 
of if 23 . 

The constancy of a n in the sigmoidal plots in Figure 4 
and the good agreement observed between the observed 
onset times and Eq. (20) indicate that the spin dynam- 
ics are dominated by the nearest-neighbor interactions in 
this regime. This is not surprising, as we are still operat- 
ing in the short time regime. Higher-order spin processes, 
if significant, would be expected to manifest themselves 
at later times, leading to a deviation from the simple 
model behavior described above. 

In conclusion, we have presented a new experimental 
method to characterize multi-spin dynamics in solid-state 
NMR free induction decay. The initial creation of coher- 
ences were observed to follow a sigmoidal growth curve, 
with the onset times characterizing the dynamics of the 
spin system. These dynamics in turn were critically de- 
pendent on the geometrical arrangement of the spins as 
expected. 

The authors thank D. Greenbaum, P. Capperallo and 
T. S. Mahesh for helpful discussions, and the NSF, ARO 
and DARPA DSO for financial support. 
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Examining Eq. (18), we see that during the initial growth 
of gK is dominated by W t K _ 1 as gx+i = 0. Thus the incre- 
mental time to go from a (K — l)-spin system to a 7f-spin 
system is inversely proportional to Wk-i- Applying this 
recursively, we obtain a general expression for the onset 
time tK = r lH 1 K ^( d , where r) is a constant and is a 
generalized Harmonic number of order r. The fits obtained 
with this equation are also in good agreement with the ex- 
perimental data, and appear to be functionally identical to 
the Eq. (20) in the range of large K. 



